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Abstract 

Milz and Strunz recently reported substantial evidence to further support the previously con¬ 
jectured separability probability of for two-qubit systems (p) endowed with Hilbert-Schmidt 
measure. Additionally, they found that along the radius (r) of the Bloch ball representing either of 
the two single-qubit subsystems, this value appeared constant (but jumping to unity at the locus 
of the pure states, r = 1). Further, they also observed (personal communication) such separability 
probability r-invariance, when using, more broadly, random induced measure {K = 3,4,5,...), 
with AT = 4 corresponding to the (symmetric) Hilbert-Schmidt case. Among the findings here is 
that this invariance is maintained even after splitting the separability probabilities into those parts 
arising from the determinantal inequality > \p\ and those from \p\ > \p^'^\ > 0, where the 

partial transpose is indicated. The nine-dimensional set of generic two-re[al]bit states endowed with 
Hilbert-Schmidt measure is also examined, with similar r-invariance conclusions. Contrastingly, 
two-qubit separability probabilities based on the Bures (minimal monotone) measure diminish with 
r. Moreover, we study the forms that the separability probabilities take as joint (bivariate) func¬ 
tions of the radii (r^, vb) of the Bloch balls of both single-qubit subsystems. Here, a form of Bloch 
radii repulsion for separable two-qubit systems emerges in all our several analyses. Separability 
probabilities tend to be smaller when the lengths of the two radii are closer. In Appendix A, we re¬ 
port certain companion analytic results for the much-investigated, more amenable (7-dimensional) 
A-states model. 
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I. INTRODUCTION 


A considerable body of diverse evidence-tliough yet no formal proof-lias been ad¬ 
duced strongly indicating that the probability that a generic two-qubit system is separa¬ 
ble/unentangled is [IHB] [71 sec. VII] [HI sec. 4]. The probability is computed with 
respect to the Hilbert-Schmidt (flat/Euclidean) measure [9l [10] on the 15-dimensional con¬ 
vex set of 4 X 4 density matrices (p). Milz and Strunz have recently conducted an analysis 
further supportive of this conjecture, while injecting an interesting new element m- They 
found that the probability of ^ appears to hold constant in the radial direction (r^) of the 
Bloch ball/sphere parameterization (r^ G [0,1], 9a G [0,27r), (j)A G [0,7r]) 


Pa = TtbP = 


( 1 ) 


1 1 cos {(j)A) r-A + l cos {9a) sin {(j)A) ta - i sin {9a) sin {(j)A) ta 

2 cos {9a) sin (0^) ta + i sin {9a) sin (0^) ta 1 - cos (0 a) ta 

of either of the qubit subsystems (A, B) of p, obtained by the partial tracing over p of the 
complementary subsystem-with a singularity occurring at the pure state boundary, ta = 1 


(cf. Fig. 10 below). (At times below, we use the symbol r, to denote interchangeably, ta 
or tb- “The Bloch sphere provides a simple representation for the state space of the most 
primitive quantum unit-the qubit-resulting in geometric intuitions that are invaluable in 
countless fundamental information-processing scenarios” [T2].) 

This same r-invariance phenomenon appeared to hold, in general they found, for 2 x m 
[qubit-qudit] systems. (For m > 3, the probability of having a positive partial transpose- 
which is now a necessary but not sufficient condition for separability-was employed m 
Fig. 5].) Further, Milz indicated in a personal communication that for the 2x2 qubit- 
qubit systems endowed with random induced measure (a function of the dimension K of 
the ancillary space) |T2l HI], r-invariance of these separability probabilities-the values of 
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which can now be directly obtained from eqnations (3)-(5) in [T5]-also seems to hold. The 
(symmetric) case = 4 is eqnivalent to the Hilbert-Schmidt one. 

The work of Milz and Strnnz is rather similar in motivation with earlier efforts in which 
it was songht to describe separability probabilities also as fnnctions of single variables (bnt 
other than the Bloch radins)-namely, the “cross-prodnct ratio” (snggested by work of Bloore 
US]). SiS [HIITIITH], and the maximal concurrence over spectral orbits [19] (cf. [20|), and 


the participation ratio and von Neumann-Renyi entropies |2ll Figs. 2b, 4] (cf. App. B below 
and [51 Figs. 1, 2]). None of the univariate separability functions constructed was of the 
highly intriguing constant form, however. (So, in a sense here, the Bloch radius serves as the 
extreme antithesis of an entanglement indicator, being entanglement-insensitive apparently.) 

In this study, we seek to broaden the investigation of Milz and Strnnz by examining the 
nature of the joint (bivariate) distribution of the Hilbert-Schmidt separability probability 

). Further, we 


over the radii {rAiTB) of both qubit subsystems of p (sec. Ill) (cf. 
examine the use of random induced measure for the Hilbert-Schmidt {K = 4)-“neighboring” 


cases of iF = 3 and 5 (sec. IV). We will, similarly, examine the separability probability with 
respect to the another measure of substantial interest, the Bures [THl EH ES] (sec. 0 , con¬ 
cluding that-contrastingly-the Bures separability probability is not constant over r, but di¬ 
minishes with the Bloch radius. Two-re[al]bit [26lEZ| Hilbert-Schmidt analyses are included 
in sec 


VT and an apparent related “Dyson-index” effect [28] in sec. VIA In sec. VH 


we 


examine how the various separability probabilities change as a function of r when they are 
subdivided into that part arising from the determinantal inequality \p^'^\ > |p| and that from 
IpI > \p^'^\ > 0 (cf. [HEH]). In that setting too, r-invariance appears to hold. Further, we 
attempted companion analyses for two-quater[nionic]bit systems [HIHEZI, but encountered 
certain conceptual/computational issues we have yet to successfully resolve. 

Appendix A develops upon the V-states analyses of Milz and Strnnz m Apps. A, B], 
finding similar results to those reported in the main body of this paper, but now also deriving 
certain exact formulas. In Appendix B we report some results (not directly pertaining to 
Bloch radii) based on an entanglement measure of Holik and Plastino [301 ^q. (9)]. 
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II. BACKGROUND 


Milz and Strunz found numerically-based evidence that the Hilbert-Schmidt volumes of 
the two-qubit systems and of their separable subsystems were both proportional to (1 — r\)^ 
[m eqs. (23), (30),(31)], with the consequent constant reitio of the two (simply proportional) 
volume functions being the aforementioned separability probability of 

For the recently much-investigated “toy” model of X-states [3TVI33] |M[ App. B], occu¬ 
pying a seven-dimensional subspace of the full hfteen-dimensional space, it was possible for 
them to formally demonstrate that the counterpart volume functions, somewhat similarly, 
were both again proportional, but now to (1—r^)^ (the square root of the higher-dimensional 
result). The corresponding (constant, but at = 1) separability probability was greater 
than A, that is | [HI Apps. A, Bj. This | result was also subsequently proven in [15], along 
with companion X-states hndings for the broader class of random induced measures [isnnj. 

However, the joint distributions over the two radii in which we are expressly interested 
here in discerning (either in the X-states and/or full model), did not seem readily derivable, 
even in the analytical frameworks of those two X-states studies (cf. [3S|)- So recourse 
to numerical methods seemed indicated. The two marginal univariate distributions of the 
desired joint bivariate volume distributions should, of course, be proportional to (1 — r\)^ 
and (1 — r\Y the full and X-states models, respectively. (In Appendix A, we do succeed 
in constructing the desired bivariate X-states total and separable volume functions and, 
hence, the associated separability probability function.) 


III. HILBERT-SCHMIDT ANALYSIS 

We generated 2,548,000,000 two-qubit density matrices, randomly with respect to Hilbert- 
Schmidt measure, using the simple (Ginibre ensemble) algorithm outlined in [361 eq. (1)]. 
(That is: “(a) prepare a square complex random matrix A of size N pertaining to the Ginibre 
ensemble (with real and imaginary parts of each element being independent normal random 
variables); (b) compute the matrix H = AA?/(tr AA?), (generically positive dehnite)” [T3].) 
For each such matrix, we found the values of r^ and r^, as well as performed the well- 
known Peres-Horodecki (determinantal-based [37]) test for separability [331ES] on the partial 
transpose of p. (That is, a two-qubit state is separable if and only if this determinant is 
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nonnegative.) Then, we discretized/binned the values of the r^’s and r^’s obtained to lie 
in intervals of length Thus, we obtained two 100 x 100 matrices of counts (which we 
symmetrize for added stability). In Figs. [^and|^ we show the histograms of these two sets 
of counts (cf. [m Fig. 3]). (Of the 10,000 bins, 9,364 of the total and 9,199 of the separable 
ones are occupied.) 

The hrst {total counts) plot appears to be somewhat broader in nature than the second 
{separable counts) plot, while appearing qualitatively rather similar. 

A natural null (product/independence) hypothesis to adopt to explain the nature of 
Figs. [^and|^in light of the assertions of Milz and Strunz HH -is that both these surfaces 
are proportional (taking into account the spherical area formula) to 

Q(r^, m) = IOttV^ r|(l - - r|)^ (2) 

As (visual) tests of these hypotheses, we show the residuals (Figs. and |^-that is, the 
discrepancies (the yet unexplained amount) from the hrst two hgures-based on such predic¬ 
tions. Q{rA,TB) was normalized so as to minimize the sum of squares of the residuals. We 
note, interestingly, that both sets of residuals are bimodal in nature-as opposed to seemingly 
random in nature-so a null product/independence hypothesis does not seem suitable, and 
further analysis is required. (We can somewhat improve the hts [more so, in the total count 
case] by employing instead Q{rA,rB) = 167r^r^r| ^ .) 

In Fig. we show the estimated joint separability probabilities (the ratio of the surface 
in Fig. 1^ to that in Fig. -which are clearly now, in contrast to the univariate case-nof 
uniform over their (unit square) domain of dehnition. The initial motivation for this study 
was to discern the functional nature of this derived (separability probability) surface. 

In Fig. 1^ for ease of visualization purposes, we perform a | rotation of Fig. g(cf. na 
eq. (7)]). 

In Fig. we present the cross-section of Fig. and indicate a closely-htting 

model for it. (Let us observe-certainly quite consistently with this hgure-that the four Bell 
states are themselves unpolarized, that is ta = tb = 0, and maximally entangled (cf. |TT]L1 
For this curve, the total volume-forming the denominator of the separability probability 
curve-appears to be proportional to (1 — r)^(l -|- 8r), and the numerator comprised of the 
separable volume to be proportional (1 — r)®(l -|- yr -|- 29r^), leading to a factor of (1 — r) 
in the equation of the curve itself. The sample separability probability for those states for 
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which ta = t_b, was recorded as 0.22753675. 

Quite contrastingly, in Fig. we show the = 1 (U-shaped) cross-section of the 

two-dimensional separability probability surface (Fig. |^. A joint plot of these last two 
[ta = tb and rA + tb = 1) curves is given in Fig. |^the first indication we have here of the 
titular “Bloch radii repulsion” effect. 

The estimated separability (marginal) probabilities over either one of the Bloch radii are 


shown in Fig. 10, along with a 95% confidence interval (based on the suitable large-sample 


normal distribution approximation to the binomial distribution) about the conjectured value 
of ~ 0.242424. (The sample estimate of this probability that we obtained here was 
0.242425003.) This plot helps to confirm the chief (r-invariance) 2x2 findings of Milz and 
Strunz [m Fig. 4]. 

IV. RANDOM INDUCED MEASURE ANALYSES 

We explore the questions raised above, but now in the broader context of random induced 
measure [ini[i3i[ii], involving the use of the natural, rotationally invariant measure on the 
set of all pure states of a 4 x 77 composite system (with K = A yielding the Hilbert-Schmidt 
measure). By tracing over the i7-dimensional ancillary system, one obtains the two-qubit 
states that we will analyze. We generate random matrices with respect to these measures 
using the algorithm specified in |32] (cf. [43] 1. 

A. The case K=3 


Setting 7 = 77 — 4 = —1, in equation (2) of the recent study 


pqubit _ 


3 4^+3(27(7 + 7) + 25)F {k + |) F(27 + 9) 


(3) 


^/7^T{3k + 13) ’ 

we obtain that the associated separability probability for this scenario is ~ 0.0714285. 
The related figures-based on 764,000,000 randomly generated 4x4 density matrices-are 
Figs. 11 to 17 (Fig. [11] total counts histogram; [12] separable counts histogram; [13] joint 


separability probability estimates; [14] = tb curve; [15] -l- = 1 curve; [16] joint 

plot of rA = tb and = 1 curves; and [17] separability probability estimates over 

Bloch radius). The sample separability probability estimate is 0.0714333, quite close to the 
theoretically-predicted value. 
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For the rA = rB curve (Fig. 14), the total volume-forming the denominator of the separa¬ 
bility probability curve-appears to be proportional to (1 —r)®, and the numerator comprised 
of the separable volume to contain a factor (1 — r)®, leading-again, as in the K = 4 case-to 
a factor of (1 — r) in the equation of the curve. The sample separability probability for those 


states for which ta = tb was recorded as 


3396373 

53960055 


0.0629424. 


B. The case K=5 


Now, inserting k = 1 into the two-qubit random-induced formula ([^ above, we obtain the 


associated separability probability for this scenario, ^ ~ 0.426573. The related hgures 


143 


parallel to those (Figs. [ll]|l7) ) for the K = 3 analysis-are now Figs. 18 to 24 (Fig. [18] 
total counts histogram; [19] separable counts histogram; [20] joint separability probability 
estimates; [21] curve; [22] ta + Tb = 1 curve; [23] joint plot of and 

fA + TB = 1 curves; and [24] separability probability estimates over Bloch radius). These are 
based on 1,267,000,000 randomly generated density matrices. The estimated separability 
probability is 0.426549. 


For the ta = tb curve (Fig. 21), the total volume-forming the denominator of the sepa¬ 


rability probability curve-appears to be proportional to (1 — r)^^(l -|- Hr -|- 40r^), and the 
numerator comprised of the separable volume to contain a factor (1 — r)^^, leading, again, 
as with 77 = 3,4, to an apparent factor of (1 — r) in the equation of the curve. 

V. BURES ANALYSIS 


Our analyses now are based on 424,000,000 randomly generated 4x4 density matrices 
with respect to the Bures (minimal monotone) measure [TOl l24] , using the algorithm given 
in [36l eq. (4)]. Paralleling those hgures given above for the K = 3 and K = 5 random 


induced measures, we present a corresponding Bures-based series of hgures (Figs. 25 -31) 


(Fig. [25] total counts histogram; [26] separable counts histogram; [27] joint separability 
probability estimates; [28] ta = tb curve; [29] = 1 curve; [30] joint plot of ta = tb 

and -|- rs = 1 curves; and [31] separability probability estimates over Bloch radius). 

Fig. [^indicates that the separability probability in the radial direction of either reduced 
qubit subsystem is not constant, but diminishes with r, in strong contrast to the cases 










analyzed above. The estimate of the Bures separability probability [HI El sg itself that 
we obtain is 0.0733096. A “silver mean” conjecture for this last value, that is ~ 

0.0733389, where the silver mean is dehned as aAg = — 1, had been advanced in |3Hl eq. 

(16)]. Clearly, however, the supporting case for this decade-old Bures two-qubit separability 
probability conjecture (based on quasi-Monte Carlo sampling, and not the more recent 
algorithm [36]) is not nearly as strong as is the multifaceted case that has been accumulating 
in the past few years for the corresponding Hilbert-Schmidt conjecture of CHg. 


VI. TWO-REBIT HILBERT-SCHMIDT ANALYSIS 


We have been noting that a remarkably strong, diverse body of evidence [IHH] -though 
yet no formal proof (cf. [TH|)-has been accumulating in the past several years for the 
proposition that the Hilbert-Schmidt separability probability of generic (15-dimensional) 
two-qubit states is ~ 0.242424. Accompanying these results has also been evidence that 
the Hilbert-Schmidt separability probability of generic (9-dimensional) two-re[al]bit states 
[261127I is If ^ 0.453125. 


In Figs. 32 38 we present a parallel set of hgures to those above, now for 4x4 density 
matrices with real entries with respect to Hilbert-Schmidt measure (Fig. [32] total counts 
histogram; [33] separable counts histogram; [34] joint separability probability estimates; [35] 
curve; [36] = 1 curve; [37] joint plot of ta = tb and = 1 curves; 

and [38] separability probability estimates over Bloch radius). (We follow the prescription 
given in [361 P- 7] regarding the generation of such random matrices, of which we generate 
2,751,000,000.) The separability probability estimate that we obtain is 0.453115. In this 
case, each Bloch radius has a two-dimensional, rather than a three-dimensional character 
(nor one-dimensional aspect, as will be the case with the X-states). 

The two 100 X 100 matrices of total counts in the Hilbert-Schmidt two-rebit analysis 
(Fig.|^ and in the Hilbert-Schmidt two-qubit analysis (Fig.[^ have 8,325 of the 10,000 cells 
both of size at least 100. The correlation between the estimated separability probabilities for 
those two sets of 8,325 cells was 0.987954. On the other hand, if we hrst square the values of 
the two-rebit separability probabilities-as random matrix (Dyson-index) theory [2S| might 
suggest-the correlation is slightly higher, 0.989462. 
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A. Dyson-index relations between rebit and qubit = rs formulas 


Further pursuing this Dyson-index ansatz, we plot in Fig. 39, the ratio of the estimated 
two-qubit = tb probability distribution to the square of its two-rebit counterpart. We 
see that this ratio appears to hold quite constant at roughly |. This ratio is also, clearly, 
formally a ratio (R = of two auxiliary quantities. The numerator Rgep is itself the 

ratio of the two-qubit separable volume to the square of the two-rebit separable volume, and 
the denominator Rtot, likewise in terms of total volumes. Our well-htting estimates of Rgep 
were 0.0023976r -I- 0.019464 and of Rtot, 0.0023868r-|-0.0171397, respectively. So, the slopes 
of the two lines are extremely close, yielding-ma the equivalence -^^-the near flatness of 
Fig.|3^ 

Along such “Dyson-index” lines of thought, we have been investigating-with some nu¬ 
merical encouragement-the possibility that in the ta = tb scenarios in the two-rebit case, 
the volume functions might have factors of (1 — r)'^ and (1 — r)5 in the total and separable 
cases, respectively. Such an investigation stemmed from the apparent occurrences, noted 
above, of factors (1 — r)® and (1 — r)® in the two-qubit counterpart volume functions-which 
appear to be, in full, proportional to (1 — r)®(l -|- 8r) and (1 — r)®(l -|- yr -|- 29r^). 

We would anticipate that related Dyson-index patterns would continue to hold in the 
case of two-quater[nionic]bits. 


VII. DIVISION OF SEPARABILITY PROBABILITIES BETWEEN > \p\ 

AND \p\ > \pP'^\ > 0 


It now appears, based on the above analyses, that for the two-qubit states endowed 
with random induced measure, the separability probabilities are constant along either Bloch 


radius (except at the isolated point r = 1) of the reduced single-qubit states (Figs. 10 


mull- An interesting supplementary question which we will now investigate is what are 
the contributions to the separability probabilities arising from the determinantal inequality 
> IpI and, complementarily, from |p| > > 0. 
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A. Hilbert-Schmidt {K = 4) case 


We know from preceding work [^l Table IV] [22] that in the Hilbert-Schmidt case, the 
conjectured probability of appears to be evenly divided, with each inequality contributing 
T. In a further analysis conducted here, based on 1,419,000,000 random matrices, it appears 
that this amount of T ~ 0.121212 (the sample estimate being 0.121208) is-paralleling the 
Milz-Strunz hnding for the total (undivided, that is > 0) separability probability- 


constant along the Bloch radius of either reduced single-qubit system (Fig. 40). 


B. Random Induced K = 3 case 

For the K = 3 {k = —1) case, the entire separability probability of T is associated with 
the inequality > |p|, so no similar nontrivial splitting can take place. (In this scenario, 
p must possess at least one zero eigenvalue, and hence |p| = 0, thus explaining this unique 
phenomenon.) 


Random Induced K = 5 case 


For the K = 5 {k = 1) case, the total separability probability appears-as already noted 


(sec. IVB)-to be 0.426573. Table IV of 0 asserts that the proportion of this associated 


with > IpI is yielding then (^)(^) = ^ ~ 0.157343. In Fig. 


41 


we show-based 


on 1,267,000,000 random matrices-an associated flat-like plot, conforming closely to this 
value (the overall sample estimate being 0.157323). Thus, it appears that this separability 
(sub-)probability-and, of course, its complementary value of ^ ~ ^ constant 

along either Bloch radius. So, these analyses serve as an expansion-and a type of further 
validation-of the Milz-Strunz hndings m- 

D. Two-rebit case 


In Fig. I^we present the Hilbert-Schmidt two-rebit counterpart of these several hgures, 
with the sample estimate of the overall separability probability of |(||) = ^ ~ 0.226563, 
being 0.226554. (Here, the 2,751,000,000 random density matrices employed in sec. VI were 
again employed.) 
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E. Bivariate Extension 


Additionally, if we similarly split the three bivariate separability probability plots (Figs. 


20, 34) for the K = A, K = 5 two-qubit and K = 4 two-rebit cases, in accordance with the 
two determinantal inequalities, \p^'^\ > |p| and |p| > > 0 , the resultant plots appear 

alike in shape to the parent plots. So, we can certainly conjecture a similar equal splitting 
of probability phenomenon in that higher-dimensional domain. 

VIII. CONCLUDING REMARKS 


All the bivariate separability probability estimates presented (Figs. ig g 1 ^ g 
appear to have a saddle point at ( 5 , 5 ), or somewhere in the neighborhood thereof, with the 


curves 


(Figs. 1^ 1^ 1^, possibly achieving their maxima at 


TA = TB = n 


and the -f- = 1 curves (Figs. ^2 29 36) attaining their minima there. A simple 

probability distribution over [ 0 , 1 ]^ with such a saddle point property, that, in addition, has 
the required marginal univariate uniform distributions over or r^, is 


pirA, tb) = 2rA + 2rs - 


(4) 


(This functional form was suggested by Brian Tung in response to a Math Stack Ex¬ 
change question https://math.stackexchange.com/questions/1271549/bivariate-probabihty- 


distributions-over-unit-square-uniform-marginals-midpo.) In Fig. 43 we show the residu¬ 


als/discrepancies obtained by subtracting from the estimated two-qubit Hilbert- 

Schmidt separability probabilities of Fig. 

We were able to obtain a somewhat superior £t (also satisfying the marginal constraints) 
to this one-as measured by the sum of absolute values of residuals from the £t-using a 
higher-degree form of probability distribution over the unit square, namely 

5rB 


!( N 3 2 ,32. Zr\ hr A 

P [Ta, Tb) = -jrA^B + i^TArB - QrArB -^ ^ 


3r| 


+ 


(5) 


4 ' 2 

From Figs. [9 16, 23, 30 and 37 we see a form of Bloch radii repulsion. That is, separability 
probabilities tend to increase as the gap in value between the lengths of the two radii increase. 

At this point, we have not yet achieved our motivating goal in undertaking this study, 
that is, to determine the precise nature of the bivariate distributions over the pair of Bloch 
radii. 


12 



























A remaining related case that is still not successfully analyzed is that of the 27- 
dimensional set of generic two-quat[ernionic]bits [27j, for which the Hilbert-Schmidt separa¬ 
bility probability appears to be 0.0804954 n An interesting question here is how 
to determine the corresponding “Bloch radii” for randomly generated two-quaterbit states 
(cf. |16]). Further, we have not yet developed a computationally feasible (Mathematica- 
implemented) algorithm for the random generation of such matrices (cf. [281 Fig. 1] [361142] ). 

Let us note that the two reduced qubit systems of a pure two-qubit system must have 
their Bloch radii equal (that is, totally “non-repulsive”). The separable pure two-qubit sys¬ 
tems form a four-dimensional submanifold of the six-dimensional manifold of pure two-qubit 
systems uni p. 368], and thus are of relative measure/probability zero. (These observations, 
it would seem, at least in an informal qualitative manner, are not inconsistent with our 
general set of results.) In terms of the “pseudo-pure” two-qubit states, that is those having 
only two distinct nondegenerate eigenvalues, Scutaru has shown that ’’the Bloch vectors of 
the corresponding qubits are related by a rotation” [47] . 

In Table we present the results of an auxiliary set of analyses. Five million random 
density matrices were generated for each scenario indicated, and the correlation computed 
between the lengths of the corresponding Bloch radii, both for all the density matrices 
generated, and also just for the subset of separable density matrices. The consistently 
smaller correlations for the separable states are a manifestation of the repulsion effect we have 
documented in this study. (We note, however, that none of these correlations is negative. 
So, perhaps rather than the term “repulsion”, the use of “relative repulsion” or “diminished 
attraction” might be more strictly appropriate.) Obviously, the correlations in the table 
based on the Bures measure are exceptionally large. 
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scenario 

all states 

separable states 

Random Induced K = 5 

0.145496 

0.0968024 

Hilbert-Schmidt two-qubits 

0.183026 

0.107762 

Hilbert-Schmidt two-rebits 

0.176898 

0.118049 

Random Induced K = 3 

0.248993 

0.125835 

Bures two-qubits 

0.388250 

0.210838 


TABLE I: Correlations between pairs of Bloch radii for all states and all separable states 
for differing scenarios. The lower correlations for the separable states are consistent with a 

relative repnlsion effect. 


Evidence addnced by Milz and Strnnz indicated that both the Hilbert-Schmidt total 
and separable nnivariate volnme fnnctions of two-qnbit states were simply proportional to 
(1 — [H], eq. (23)]. Qnite snpportingly, we ourselves £t a function of the form c(l — r'^y 


to the large sample of such states employed above (sec. Ill), and obtained estimates of 
5.99965 and 5.99926, respectively, of this exponent for the two volumes. Similarly, for the 


two-rebit set of analyses (sec. VI), the estimates were 6.00439 and 6.00447. For random- 


induced measure with K = 3 (sec. IVA), 3.99973 and 4.00015 were obtained, while for the 


case K = 5 (sec. IVB), the corresponding results were 7.99923 and 7.99917. (A parallel 
exercise based on the Bures measure [sec. |V] yielded the rather proximate-but certainly far 
from integral-results of 3.48845 and 3.58319, respectively.) 

Following the work of Braga, Souza and Mizrahi [iQl eq. (7)], it might prove advantageous 
in our quest to model the various bivariate total and separable volume and probability 
functions discussed above, to employ transformed variables of the form and 

u_ = In fact, this appears to be the case in the following appendix devoted to 

V-states analyses 
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IX. APPENDIX A: A-STATES ANALYSES 


A-states bivariate formulas 


We employed the A-states parametrization and transformations indicated by Braga, 
Souza and Mizrahi HQl eqs. (6), (7)]. Based on these, we were formally able to repro¬ 
duce the Hilbert-Schmidt univariate volume result of Milz and Strunz [m eq. (20), Fig. 


Vhs - 2304* 




( 6 ) 


(7) 


as the marginal distribution (over either or r^) of the bivariate distribution (Fig. 44) 

{^a - 1) ^ {ta {ta + 3) - 5r| -hi) > tb 
- 1) ^ (-5ri -h Tb (vb -h 3) -h 1) ta < vb- 

To now obtain the desired A-states bivariate separability probability distribution (perhaps- 
and hopefully-suggestive of the full 15-dimensional counterpart), we hnd the separable vol¬ 
ume counterpart (Fig. 45) to ([^ 

sepV^s\rA,rB) = 

’r^(»’A-l)^(5(r'A-l-3)r|-10(3rA+l)r'|+8r^+9rA+3) 


( 8 ) 


7680 


Ta > Tb 

7r^(»’s-l)®(5r^(rs+3)-10r3^(3rs+l)+rs(8rB+9)+3) ^ ^ ^ 

7680 < ^B, 


and take their ratio (Fig. 46) (noting the cancellation of the (r — l)^-type factors), obtaining 
thereby the A-states bivariate separability probability formula. 


5(rA+3)r|j-10(3rA+l)r-|+8r^+9rA+3 
8(r’A(rA+3)-5r|+l) 
5r\{rB+3)-l0r\(3rB+l)+rBi8rB+^)+3 


(9) 


Ta > Tb 
Ta < Tb 


8 (-5r^+rs (rs+3)+1) 

(Numerical integration of this function over [0,1]^ yielded 0.381678.) Also, Fig. 47 shows 


the (lower) ta = tb and (upper) ta + tb = 1 cross-sections of Fig. 46, (We note in 


the next section immediately below an exception to the “repulsion” phenomenon, which 
we have repeatedly observed above. We computed the correlation between the pair of 
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Bloch radii to be 1 — 3774 S 736 _ioogo^ 2 _|_^4 ~ 0.702341 for all states and only slightly less, 
1 “ 235929600 ^- 252007 r 2 + 7 r 4 ~ 0.68326, for the Separable X-states.) 

In light of the X-states results ([^ and ([^, we might speculate that the counterpart 
bivariate total and separable volumes for the 15-dimensional set of two-qubits states will both 
consist of the product of (1 — r)® and certain polynomials. The corresponding separability 
probability function (cf. Fig. would then be a rational one. 


B. Certain univariate X-states separability probability conditional distribntions 


The analytic form of the ta = tb separability probability curve for the X-states is 


p{X-states}i^^^ = Tb) 


(r - l)(5r(r(r+ 5) + 3) + 3) 
32r -I- 8 


( 10 ) 


The value of this X-states separability probability univariate function at (|, |) is ||| = 
at (0, 0) it is |, and at (1,1) it is 0. The maximum of the ta = tb curve is achieved at the 
positive root (r 0.2722700792) of the cubic equation 3r^ -|- 9r^ -|- r — 1 = 0, its value there 
(~ 0.393558399) being the positive root of the cubic equation 54r^ -|- 108r^ — 28r — 9 = 0. 
On the other hand, the minimum (||| 0.361979) of the rA + rB = l curve 

{ (rA-2)rA(5rA(?-^+r-A-10)+28)-|-8 
8(rA(4r-A-13)+4) 

rA(r-A(5rA((rA-4)rA-6)+32)+25)-20 . 

8(rA(4rA+5)-5) ^ 4 


is attained more simply at = 4. The maximum of 4 is situated at = 0 or 1. 

So, at least in this model there does not seem to be a corresponding minimax result. (Let us 
very interestingly note that in the interval r G [0.40182804, 4], the p{x-states} 
curve is dominated by the p{x - states} The maximum gap of 0.0056796160 

is attained at r = 0.4564893379. So, at least in the X-states setting, the ’’repulsion phe¬ 
nomenon” does not fully hold.) 

Further, setting = 4^ we have 


p{X - states} 


Putting = 0, we obtain 

p{X - states} 


rA(128rA+29)+23 1 ^ ^ i 

32(4rA(rA+3)-l) 2 ^ ^ ^ ^ 


35rjj-50r^+19 

44-80r^ 


0 < Ta < I 


_ Ta (8rA 9) -F 3 
^“8(rA(rA + 3) + l)’ 


( 12 ) 


(13) 
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and with tb = 1, 


= 1) = 0. (14) 

It certainly appears that this last resnlt is in direct contradiction with certain assertions of 
Milz and Strnnz: “The latter fact that psep {l) = 1 is clear: a pnre rednced state (r = 1) 
can only be realized by a prodnct and thus, a separable total state” m p. 8] (see also the 
discussion prior to their eq. (23)). We anticipate an eventual clarihcation of this apparent 
conflict, possibly in terms of differing dimensionalities of the measures employed. 

If we restrict the X-states to those for which = tb, then the separability probability 
for this continuum r/^ G [0,1] of states is ^ ~ 0.3809524. For those for which = 1 — r^, 
the corresponding separability probability is slightly higher, ^ 0.3945578. 


C. Use of Fano correlation parameter C 33 = 


In their X-state studies, both Milz and Strnnz mi and Braga, Souza and Mizrahi 001 
employ the well-known Fano parameterization of two-qubit systems [IH]. Milz and Strnnz 
denote the Fano correlation parameter in the (conventionally denoted) or a; 3 -direction 
by C33, while Braga, Souza and Mizrahi employ the notation M^z- (The alignments of the 
Bloch radii-r^ and r^-are along this same direction in the X-states model, we interestingly 
note.) Focusing on this parameter yields a number of analytic results, such as the associated 


X-states separability probability (Fig. 48). (Numerical integration of this function over [-1,1] 
yielded 0.416283.) 


D. Random Induced K = 5 case 


We have found here-introducing a factor of \p\ into the integrations in the previously- 
conducted Hilbert-Schmidt X-states analyses-that the total volume bivariate distribution 
for the induced measure case of X = 5 equals 


totVl^tl{rA,rB) = 

{ 7r2(rA-l)®(-6rA(m+5)r|+rA(rA+l)(m(m+4)+5)+21r|j-6r|+l) 

1290240 

7r^(»’s-l)®(-6r-^()’s(rs+5)+l)+21r^-|-rfl(rs+l)(rs(rs+4)+5)-|-l) 

1290240 


Ta > Tb 
Ta < Tb 


( 15 ) 
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The total volume itself is ggjg 2 oo ~ 9.8901759 ■ 10 The marginal distributions of the total 
volume bivariate distribution (15) are of the (again, proportional to 7r^(l — form 

^2(1 2)5 ^2(1 2)5 


= 


3686400 


214 . 32.52 


The separable volume is given by 


sepVlS^irA.rB) = 


7r^hA-l)®(-21(rA+5)r|+63(5rA + l)r|-27r^(8rA+5)r|+r^(8rA(rA+2)(rA+3)+25)-27r|+5) 

10321920 

7r^(?’s-l)®(-21r-^(rB+5)+63r^(5rs+l)-27r^(rs(8rs+5)+l)+25rB+8r|(rs+2)(rB+3)+5) 

10321920 


(16) 

(17) 

Ta > Tb 
ta < Tb 


Its marginal distributions are of the form 


sepVx=5\' ) — 


^2(1 2)5 ^2(1 2). 


5734400 2^5 ■ 52 ■ 7 ■ 

The separability probability we found was ^ ~ 0.642857-a result also derivable from a 
formula [IHl P-13], 

Pr{|p^'^| > 0} = 1 


2T{2k + Af 


r(fc + 2)r(3fc + 6)’ 

inserting k = 1. The separability probability function (Fig. 49) is 

vl-jr^^TA^TB) = 

-21(rA+5)r|+63(5rA+l)r|j-27rA(8rA+5)r|+rA(8rA(rA+2)(rA+3)+25)-27r|+5 
8(-6rA(rA+5)r|+rA(r'A+l)()'AhA+4)+5)+21r|-6r|+l) 
-21r^(rB+5)+63r^(5rB + l)-27r^(rfl(8rfl+5)+l)+25rB+8r|(rfl+2)(rB+3)+5 
8(-6r^(rs(rs+5)+l)+21r^+rs(rs+l)(rs(rs+4)+5)+l) 


(19) 


( 20 ) 


Ta > Tb 


Ta < Tb 


In Fig. 50, we show the ta = tb and = 1 — sections of this plot. The minimum of the 


Ta = 1 — Tb section is once again found at r = ^ with a value of = 


13-97 

21-17 


there. For the ta = tb section (cf. (10)) 


^X—states ^ \ _ 

Pk =5 yA,rB){rA=rB) — 


(1 - r)(r(21r(r(r + 8) + 6) + 40) + 5) 


8(r(16r + 7) + l) 
The maximum, 0.63964, of this curve is attained at r = 0.238465. 


0.579504 


( 21 ) 


E. Random Induced K = 6 case 


Here the total volume itself is 


9081072000 


are of the form 




a 1.0868325 ■ 10 The marginal distributions 

.2)7 ^2(i_^2)7 


2890137620 2^8 ■ 3^- 5^- 7^' 


( 22 ) 


The separability probability is given by (19), inserting k = 2. 
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F. Random Induced K = 7 case 


Here the total volume itself is 


5866372512000 

are, continuing the observed pattern, of the form 


1.68240328-10 The marginal distributions 


^ - 1)^ ^ 

’ 1664719257600 224. 34 . 52.72 ' 


(23) 


The separability probability is given by (19), inserting k = 3 there. 


Random Induced K = 3 case 


The total volume is | 7 r^log^( 2 ) 3.16125412. The bivariate total volume distribution is 


totVl^%{rA,rB) = 


-277^ log^(2) (r^ — 1) A + rs > 0 A < 1 

—277^ log^(2) {tb — 1) A > 0 A < 1 


(24) 


The two marginal distributions are of the form vr^ (1 — r^) log^(2). Here we found the sep¬ 
arability probability to equal 7 . We have not been able to hnd analytic formulas for the 
bivariate separability volume function and the bivariate separability probability function. 


but in Fig. 51, we present a numerically-based estimate of the separability probability func¬ 
tion. 


X. APPENDIX B: SEPARABILITY PROBABILITIES AS A FUNCTION OF 

\\p- ® P^Whs 


In a recent paper of Holik and Plastino, the expression 


||p-p ®p Has 


(25) 


is put forth as a measure of entanglement m eq. (9)], where 


norm 


||2l||is = tr(2lV), 


\hs is the Hilbert-Schmidt 


(26) 


In Fig. we show estimates of the two-qubit separability probabilities as a function 
of this term for the Hilbert-Schmidt and Bures measures, based on 54,000,000 random 
realizations in the former case and 43,000,000 in the latter (cf. [ 6 l Fig. 2]). The probability 
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diminishes as the Hilbert-Schmidt distance from product states ® p^) increases. (Of 
course, it might be of some interest to employ Bures counterparts of (|25| and (26)). 
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FIG. 1: Histogram of Hilbert-Schmidt randomly sampled two-qubit density matrices 

parameterized by and 
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FIG. 2: Histogram of Hilbert-Schmidt randomly sampled separable two-qubit density 

matrices parameterized by ta and tb 
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FIG. 3: Residuals (discrepancies) from a fit to the total counts (Fig. 1) of a normalized 
form of Q(r^, Tb) = 167r^r^r|(l — r^)®(l — 
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FIG. 4: Residuals (discrepancies) from a fit to the separable counts (Fig. 2) of a 
normalized form of Q{rA,rB) = 16n^r\r‘^{l — r^)®(l — r^)® 
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sep. 


FIG. 5: Estimated joint Hilbert-Schmidt two-qubit separability probabilities-the ratio of 


the separable counts in Fig. 


to the total counts in 
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sep. prob. 



FIG. 7: Estimated Hilbert-Schmidt separability probabilities for = vb- A closely-fitting 
model (ratio of apparent separable and total volumes) for this curve is 4096 ^( 8 rHn)’^^'*~^^ • 
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FIG. 8: Estimated Hilbert-Schmidt two-qubit separability probabilities for = 1. A 

closely-fitting model for this curve is 24^+8 (0 < < i) 24^-32 
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FIG. 9: Joint plot of the last two figures 
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sep. prob. 
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FIG. 10: Estimated (marginal) Hilbert-Schmidt two-qubit separability probabilities along 
either Bloch radius, together with 95% conhdence limits about the conjectured value of 

^ 0.242424 
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FIG. 11: Histogram of randomly sampled (with respect to the random induced K = 3 

measure) two-qubit density matrices 
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FIG. 12: Histogram of randomly sampled (with respect to the random induced K = 3 
measure) separable two-qubit density matrices 
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0.2 0.4 0.6 0.8 1.0 

FIG. 14: Estimated random induced {K = 3) two-qubit separability probabilities for 
'^A = ^B- A closely-fitting model (ratio of apparent separable and total volumes) for this 
curve is j^(l — r) (r^ -|- 6r -|- 1), giving a total separability probability for the continuum of 
such states of Ag 0.0629371. The maximum of the curve ^(5^10 — 14) 0.0670885 

occurs at r = |(2-\/lb — 5) 0.441518 
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FIG. 15: Estimated random induced (K = 3) two-qubit separability probabilities for 

rA + rB = 1 - 
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FIG. 16: Joint plot of the last two figures 
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FIG. 17: Estimated random induced {K = 3) two-qubit separability probabilities over 
either one of the Bloch radii, along with 95% conhdence limits about the conjectured value 

of jb ^ 0.0714285 
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FIG. 18: Histogram of randomly sampled (with respect to the random induced K = 5 

measure) two-qubit density matrices 
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FIG. 19: Histogram of randomly sampled (with respect to the random induced K = 5 
measure) separable two-qubit density matrices 
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sep. 


FIG. 20: Estimated joint random induced {K = 5) separability probabilities-the ratio of 

Fig. 


19 


to Fig. 


18 


sep. prob. 



FIG. 21: Estimated random induced {K = 5) two-qubit separability probabilities for 
= 'T'b- a closely-fitting model (ratio of apparent separable and total volumes) for this 


curve is 



40r^+llrH-l 
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FIG. 22: Estimated random induced {K = 5) two-qubit separability probabilities for 

Ta + Tb = 1 - 
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FIG. 23: Joint plot of the last two figures 
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FIG. 24: Estimated random induced {K = 5) two-qubit separability probabilities over 
either one of the Bloch radii, along with 95% conhdence limits about the conjectured value 

of ^ 0.426573 
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FIG. 25: Histogram of randomly sampled (with respect to the Bures measure) two-qubit 

density matrices 
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FIG. 26; Histogram of randomly sampled (with respect to the Bures measure) separable 

two-qubit density matrices 
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FIG. 27: Estimated joint Bures two-qubit separability probabilities-the ratio of Fig. 

Fig. 


26 


to 


25 
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FIG. 28: Estimated Bures two-qubit separability probabilities for ta = tb 
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FIG. 30: Joint plot of the last two figures 


38 






sep. prob. 
0.078 - 


0.076 - • 



FIG. 31: Estimated Bures two-qubit separability probabilities over either one of the Bloch 
radii, along with ill-fitting (in contrast to Hilbert-Schmidt and random induced) 95% 
conhdence limits about the conjectured overall probability of ~ 0.0733389. 



FIG. 32: Histogram of randomly sampled two-rebit density matrices 
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FIG. 33: Histogram of randomly sampled separable two-rebit density matrices 
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FIG. 34: Estimated joint Hilbert-Schmidt two-rebit separability probabilities-the ratio of 


Fig. 33 to Fig. 


32 
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FIG. 35: Estimated Hilbert-Schmidt two-rebit separability probabilities for ta = tb 
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FIG. 36: Estimated Hilbert-Schmidt two-rebit separability probabilities for = 1 

sep. prob 



FIG. 37: Joint plot of the last two figures 
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FIG. 38: Estimated (marginal) Hilbert-Schmidt two-rebit separability probabilities over 
either one of the Bloch radii, along with 95% conhdence limits about the conjectured value 

of i 0.4531250 
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FIG. 39: Ratio of the separability probability square of the 

separability probability = tb) 
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FIG. 40: Estimated two-qubit separability probability for the Hilbert-Schmidt case {k = 0) 
associated with the determinantal inequality |p^^| > |p|, along with 95% confidence limits 

about ^ ~ 0.121212. 
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FIG. 41: Two-qubit separability probability for the random induced measure case of 
K = b {k = 1) associated with the inequality > |p|, along with 95% conhdence limits 
about ^ ~ 0.157343, with the complementary probability for |p| > |p^^| > 0 being 


FI _= j_ ~ n 9609'?! 

143 286 26 O.zuazoi. 
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FIG. 42: Two-rebit separability probability for the Hilbert-Schmidt case associated with 
the inequality \p^'^\ > |p|, along with 95% confidence limits about 
(i)(i) = fl«0'2265625. 
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FIG. 43: Residuals (discrepancies) of £t of ^^^^^^-given by (4)-to the estimated 
Hilbert-Schmidt two-qubit separability probabilities given in Fig. 
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FIG. 45: Bivariate Hilbert-Schmidt separable volume distribution (8) for the X-state model 
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FIG. 46: Bivariate Hilbert-Schmidt separability probability distribution (9) for the X-state 


model-the ratio of Fig. 45 to Fig. 


44 


X-states sep. prob. 



FIG. 47: (Lower) ta = rs-given by (lO)-and (upper) = 1 curves-given by (ll)-for 


bivariate Hilbert-Schmidt X-states separability probability distribution. The minimum of 
the upper curve is at r = while the maximum of the lower curve is at 0.27227007. In the 
interval r G [0.40182804, 4] the ( 7 -^ = i _ 7 -^) curve is dominated by the 


P 


{X—states} / _ 


{ta = tb) curve. 
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FIG. 48: X-states separability probability as a function of the Fano correlation parameter 
C 33 = M^z- For —1 < Mzz < 0, the curve is simply 2 {m^- 2 )^^m -i) - maximum 
0.79819147) occurs near C 33 = M^z = 0.03677089. 



FIG. 49: Bivariate separability probability function (eq. (20)) for the K 

measure X-states model 


5 induced 
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FIG. 50: The ta 


tb and (largely dominant) = 1 — sections of the K 


separability probability fnnction (20) 
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FIG. 51: Nnmerically-based estimate of bivariate separability probability fnnction for the 

K = 3 indnced measure X-states model 
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FIG. 52: Two-qubit separability probabilities as a function of the entanglement expression 
\\p - ® P^Whs for the Hilbert-Schmidt (dominant curve) and Bures measures 
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